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Abstract — An extended version of the Cascade-Exciton Model (CEM) of nuclear reac- 
tions is applied to analyze nucleon-induced fission cross sections for 209 Bi and 208 Pb nuclei 
in the 45-500 MeV energy range. The available data on linear momentum transfer are an- 
alyzed as well. The results are compared with analytical approximations resulting from a 
comparative critical analysis of all available experimental data. Systematic discrepancies 
between calculations and experimental data are revealed. A modification of the CEM is 
proposed, which significantly improves the model predictions for projectile energies above 
100 MeV. 



I. INTRODUCTION 



A large amount of nuclear reaction data, including fission cross sections at intermediate 
energies is required for applications, e.g., for accelerator transmutation of waste (ATW) 
for elimination of long-lived radioactive wastes with a spallation source, accelerator-based 
conversion (ABC) aimed to complete the destruction of weapon plutonium, accelerator- 
driven energy production (ADEP) which proposes to derive fission energy from thorium 
with concurrent destruction of the long-lived waste and without the production of weapon- 
usable material, for accelerator production of tritium (APT) etc. (see the Proceedings 
of the international conferences on Accelerator-Driven Transmutation Technology held in 
Kalmar [JTp and Dubna 0] and references therein). Experiments to measure these data are 
costly and there are a limited number of facilities available to make such measurements. 
Therefore reliable models are required to provide the necessary data. 

During the last two decades, several versions of the Cascade-Exciton Model (CEM) [TJ] 
of nuclear reactions have been developed at JINR, Dubna (for an overview, see [|]]). A 
large variety of experimental data for reactions induced by nucleons [Bp, pions [^p, and 
photons ^ has been analyzed in the framework of the CEM and the general validity of this 
approach has been confirmed. The recent International Code Comparison for Intermediate 
Energy Nuclear Data [J8p has shown that the CEM adequately describes nuclear reactions 
at intermediate energies and has one of the best predictive powers for double differential 
cross sections of secondary nucleons as compared to other available models (see Tabs. 5 
and 6 in the Report [||] and Fig. 7 in Ref. [j^]). In the last few years, the CEM has been 
extended Jl^j ] to calculate hadron-induced spallation and used to study [JTTJ]— [[T3p about 700 
reactions induced by protons from 10 MeV to 5 GeV on nuclei from Carbon to Uranium. 

A detailed description of the CEM may be found in Ref. and of its extended version, 
as realized in the code CEM95, in Refs. [ |1T| Op; therefore, we mention here only its 
basic assumptions. The CEM assumes that reactions occur in three stages. The first stage 
is the intranuclear cascade in which primary and secondary particles can be rescattered 
several times prior to absorption by, or escape from the nucleus. The cascade stage of the 
interaction is described by the standard version of the Dubna intranuclear cascade model 
(ICM) [|T5p. The excited residual nucleus remaining after the emission of the cascade 
particles determines the particle-hole configuration that is the starting point for the second, 
preequilibrium stage of the reaction. The subsequent relaxation of the nuclear excitation 
is treated by an extension of the Modified Exciton model (MEM) flP7| ] of preequilibrium 
decay which also includes the description of the equilibrium evaporative third stage of the 
reaction. 

Recently, the CEM has been extended by taking into account the competition between 
particle emission and fission at the compound nucleus stage JT8| ] and a more realistic cal- 
culation of nuclear level density [[T9p. An earlier extended version of the CEM, as realized 
in the code CEM92, was used by Konshin [ |T5| ] to calculate nucleon-induced fission cross 
sections for actinides in the energy region from 100 MeV to 1 GeV. Nevertheless, pre- 
viously the CEM has not been applied to study fission cross sections for pre-actinides; 
therefore, its predictive power and applicability to evaluate arbitrary fission cross sections 
was unknown. 
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In order to investigate the applicability of the CEM to fission cross sections and hoping 
to learn more about intermediate-energy fission and to identify possible improvements to 
the CEM and other models to improve their predictive power, we use here an extended ver- 
sion of the CEM (see the recent report and references therein), as realized in the code 
CEM95 pCp, to perform a detailed analysis of proton- and neutron-induced fission cross 
sections of 209 Bi and 208 Pb nuclei and of the linear momentum transfer to the fissioning 
nuclei in the 45-500 MeV energy range. 

The 208 Pb and 209 Bi target nuclei and this energy interval were chosen for this study 
for the following reasons: 

• new experimental information on 208 Pb(n,f) and 209 Bi(n,f) reactions has been re- 



cently obtained at these incident energies (see Q2lQ-[]24|] and references therein), 



the reactions under study have practical importance in connection with concepts 
of accelerator-driven transmutation and tritium production technologies, which may 
include the use as a neutron source of massive lead or lead-bismuth targets irradiated 
by proton beams, 

cross sections for these reactions are widely used as standards in the intermediate 



energy region (e.g., p5| , [260), with cross sections for other targets being normalized 
to them, 

fissioning nuclei in this mass region are close to the doubly magic nucleus 208 Pb, 
where one would expect the clearest manifestation of shell effects, which are impor- 
tant for the development of any model, 

fission barriers of nuclei in this mass region do not have the double-humped structure 
which appears in the actinide nuclei and makes calculations for them more compli- 
cated and uncertain. 



II. CALCULATIONAL TECHNIQUE 

As the fission cross sections for nuclei under study constitute a small part of the total 
reaction cross section, the statistical weight method pTS| , |2"7[ ] is used to calculate fission 
cross sections instead of the straightforward Monte-Carlo technique. It allows us to obtain 
statistical uncertainties of the calculated fission cross sections, as a rule, not more than 
5-10% with only 3000 inelastic events in each calculation. 

In the course of the calculations, the characteristics of the fissioning nuclei (the charge, 
the mass, the excitation energy, all components of the linear and angular momenta, and the 
statistical weight of the event) are saved, and later used for the construction of distributions 
of fission events with respect to specified parameters as well as for the calculation of the 
corresponding average values. Examples of these distributions are shown in Figs. 1 and 2. 

Fig. 1 represents the distributions of calculated longitudinal linear momentum transfer 
(LMT) to the fissioning nuclei in the 209 Bi(p,f) reaction for incident proton energies of 45 
MeV (a), 73 MeV (b), 96 MeV (c), and 160 MeV (d). Comparing these distributions allows 
one to follow the changes in the reaction mechanism with increasing incident particle 
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energy. The calculation predicts that most nuclei undergoing fission induced by 45 MeV 
protons possess a linear momentum equal to that of the incident proton (Fig. la), i.e., the 
full LMT mode predominates, which corresponds to fission from the compound nucleus. 
(Hereinafter the term "compound nucleus" implies the nucleus with the charge and the 
mass equal to the sum of the corresponding characteristics of the target nucleus and of the 
incident particle.) With increasing incident proton energy, the share of the events with the 
full LMT falls quickly (Fig. lb). For proton energies of 100 MeV or more (Fig. lc, d) the 
fission from the compound nucleus practically "dies out"; the maximum in the distribution 
becomes displaced to the region of incomplete LMT and the distribution width increases. 
A comparison of the LMT calculations with experimental data is discussed in Sect. III. 

Fig. 2 shows the distributions of the charge of the fissioning nuclei (a, b), mass (c, d), 
and excitation energy (e, f) in the 209 Bi(p,f) reaction for proton energies of 45 and 160 
MeV. The curves reflect the same tendency as discussed above. At 45 MeV, fission from 
the compound nucleus ( 210 Po) predominates. On the contrary, at 160 MeV the distributions 
become displaced and their shape approaches symmetry. 

In Fig. 3, we show the dependence of the average charge (a), mass (b) and excitation 
energy (c) of the fissioning nuclei in the 209 Bi(p,f) and 209 Bi(n,f) reactions as functions 
of the incident particle energy in the 45-160 MeV range. The energy dependence of the 
average mass and excitation energy of the fissioning nuclei are nearly the same for these 
two reactions. On the other hand, the average charge of the fissioning nuclei for the (p,f) 
reaction decreases faster than for the (n,f) one. This can be qualitatively explained by 
noticing that there is a significant probability that the incident nucleon leaves the nucleus 
after an inelastic scattering during the intranuclear cascade stage of the reaction. 

The results presented in this and subsequent sections are obtained with input parame- 
ters of CEM95 as discussed in Sect. V. 
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Fig. 1. Calculated distributions of longitudinal linear momentum transferred to the fission- 
ing nuclei in the 209 Bi(p,f) reaction for incident proton energies of 45 MeV (a), 73 MeV (b), 
96 MeV (c), and 160 MeV (d). 
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Fig. 2. Calculated distributions of fissioning nuclei in the 209 Bi(p,f) reaction vs. charge (a,b), 
mass (c,d), and excitation energy (e,f) for incident proton energies of 45 MeV (a,c,e) and 160 MeV (b,d,f)- 
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Fig. 3. Calculated average charge (a), mass (b), and excitation energy (c) of fissioning nuclei 
in the 209 Bi(p,f) reaction as a function of projectile energy. 
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III. LINEAR MOMENTUM TRANSFER 



It has been previously shown that the momentum transferred in nucleon reactions 
with actinide nuclei is better described by the CEM than by other models (see details 
inRef. [|28|]). 

Fig. 4 shows the calculated average LMT for the 209 Bi(p,f) (a), 209 Bi(n,f) (b), 208 Pb(p,f) (c), 
and 208 Pb(n,f) (d) reactions in comparison with the sparse available experimental data fl29|]- 
[J32[j. A large spread exists between the 209 Bi(p,f) LMT data from different authors (Kowal- 
ski Stephan pD[], Shigaev pT|]) and even within the same experimental dataset plf]. 
This may indicate the presence of some omitted systematic uncertainties in that work []3ltl, 
which was performed with solid-state nuclear track detectors. Possible sources of system- 
atic uncertainty may include a decrease of detection efficiency at back angles, fluctuations 
of the position of the beam at the target, as well as inhomogeneity of the target. We regard 
as more reliable the data from the works [ j29] , |30| ] where fission fragments were registered 
by semiconductor detectors in coincidence mode. Judging from these data, one can con- 
clude that our present CEM95 calculation predicts correctly a decrease with energy of 
the average LMT for the 209 Bi(p,f) reaction but underestimates its average value by 1-1.5 
times the experimental uncertainty for the highest energy. 

To the best of our knowledge, only one experimental work is available on LMT for the 
(n,f) type of reactions [ f32| ] (see Fig. 4b). In this work the average LMT was obtained for 
the 209 Bi(n,f) reaction at 75 MeV by measuring the difference between the kinetic energies 
of the fission fragments registered by an ionization chamber in the forward and backward 
hemispheres. In contrast to the 209 Bi(p,f) results, the average LMT calculated with the 
CEM95 code lies above the sole experimental point by about 1.5 times the experimental 
uncertainty. 

The works of Kowalski [ |Z9"| ] and Stephan p3Dp also include measurements of LMT 
distributions for the fissioning nuclei in the 209 Bi(p,f) reaction at 156 and 96 MeV, re- 
spectively. To a first approximation, the distributions are Gaussian-shaped, similar to the 
CEM95 predictions (see Fig. lc, d). For a quantitative comparison between the calcula- 
tions and experiment, we calculate the second moment (dispersion) of the LMT experi- 
mental and calculated distributions. In order to compare the experimental dispersion with 
the calculated one, it is necessary to introduce a correction to the latter, which takes into 
account the prompt neutron emission from the fission fragments. We estimate this correc- 
tion using the average kinetic energy of the prompt fission neutrons from Ref. [ [33| ] and the 
average number of neutrons per fission from Ref. []34jJ. The comparison between calculated 
and experimental results is shown in Table I. 

Table I shows that there is good agreement between the experimental data and our 
results corrected for the prompt neutron emission from the fission fragments. 

On the whole, we conclude that the description of linear momentum transfer processes 
provided by CEM95 does not contradict the available experimental data. 

In Sect. VI we discuss a sensitivity study of the calculated fission cross sections to dif- 
ferent input parameters of CEM95. We also study the sensitivity of the LMT description. It 
is found that neither the average nor the dispersion of the LMT varies beyond the statistical 
uncertainties of the calculations for any reasonable variation of the input parameters. 
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Fig. 4. The average share of the linear momentum transferred to the fissioning nuclei as a func- 
tion of the projectile energy, for the following reactions: 209 Bi(p,f) (a), 209 Bi(n,f) (b), 208 Pb(p,f) (c), 
and 208 Pb(n,f) (d). The points are experimental data from Refs. [^]-[32]. The dashed lines and 
the open squares represent the present calculations. 
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This insensitivity can be interpreted as follows. The momentum losses occur mainly at the 
cascade stage of the reaction due to forward-peaked emission of fast cascade particles. As 
the excited nucleus nears statistical equilibrium, the secondary particles become less en- 
ergetic and their emission becomes more isotropic in the center-of-mass system. Thus the 
secondaries being emitted at the later reaction stages take away less and less momentum. 
Since the free parameters available in CEM95 relate only to these later reaction stages, it is 
natural that their variations have only a weak influence on the resulting LMT distributions. 



TABLE I 

Dispersion of the LMT distribution in the 209 Bi(p,f) reaction 
(in units of the incident proton linear momentum) 



Proton energy 


CEM95 results 


Correction due to 


CEM95 


Experimental data 


(MeV) 


(uncorrected) 


prompt neutron 


corrected 








emission by 


results 








fission fragments 






96 


0.25 


0.28 


0.38 


0.36 ±0.02 [30] 


160 


0.32 


0.21 


0.38 


0.38 ±0.02* [|290 



* The experimental value was obtained at the proton energy of 156 MeV. 



Furthermore, no notable differences are seen between the calculated LMT data for in- 
cident protons and neutrons of the same energy. The calculated results for 208 Pb and 209 Bi 
target nuclei are also similar. Thus, CEM95 predicts a weak sensitivity of the LMT pro- 
cesses to variations of the incident particle charge as well as to the target nucleus charge 
and mass. 



IV. APPROXIMATING EXPERIMENTAL FISSION CROSS SECTIONS 

In contrast to the poor experimental database for LMT, the measurements of fission 
cross sections for 208 Pb and 209 Bi are numerous although their reliability varys greatly. 
Thus the most effective way to estimate the quality of the model predictions is a compar- 
ison of the calculated results with approximations built on the basis of the experimental 
data taking into account differences in their systematic and statistical uncertainties. 

Prokofiev and co-authors [ [35] ] performed a literature search and a comparative critical 
analysis, ranking the experimental data on the (p,f) cross sections for a wide range of nuclei 
from Ta to U in the energy region up to 30 GeV. We consider only the data for 208 Pb and 
209 Bi up to 500 MeV. An approximation for the 209 Bi(n,f) cross section was suggested by 
Smirnov and published in Ref. p6|]. The approximation given below for the 208 Pb(n,f) 
cross section is obtained here and has not been published previously. 

A list of publications containing the experimental data on the nucleon-induced fission 
cross sections for 208 Pb and 209 Bi is given in Table II. 
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TABLE II 

References to publications containing experimental data on the nucleon-induced fission 
cross sections for 208 Pb and 209 Bi at energies up to 500 MeV 



Target nucleus 


Incident particle 




P 


n 


209 Bi 
208 pb 


[29, 3§, [36]-[5C] 
[©[^,[39,18,0 


[21]-[0],l]5j,gl 
H5_l,|5J,[gTl]-0 



The experimental data from the publications listed in Table II were subjected to a crit- 
ical comparative analysis [ |35| ] in order to establish how reliably the quantities defining 
fission cross sections (incident particle flux, target thickness, fission fragment detection 
efficiency, etc.) were measured. As a result of this analysis, "reliability grades" were as- 
signed to every dataset based on the following criteria: 

• the 1st (highest) grade was assigned to a dataset if the analysis did not reveal any 
omitted or underestimated sources of systematic uncertainties. Thus, the uncertain- 
ties of the lst-grade data were taken from the original publications, 

• the 2nd (medium) grade was assigned to a dataset if there were some omitted or 
underestimated sources of systematic uncertainties that could be deduced. Thus, 
enhanced uncertainties were assigned to the 2nd-grade data, 

• the 3rd (lowest) reliability grade was assigned to a dataset if the analysis revealed 
defects leading to large and unpredictable systematic errors. The 3rd-grade datasets 
were removed from further consideration. Data remeasured or reprocessed by their 
author(s) in later publications were also not considered. 

The reliability grade of a dataset was lowered if an original publication does not contain 
a detailed description of the experimental technique. 

The highest grade was assigned, as a rule, to the datasets obtained by modern tech- 
niques: ionization chambers (for the (n,f) reactions), semiconductor detectors, thin film 
breakdown counters, solid state nuclear track detectors with good geometrical conditions. 
On the other hand, the datasets obtained by the nuclear photo-emulsion technique or by 
fission product yield summation generally received the lowest grade. 

In a few cases, two or more datasets show a similar energy dependence of the cross 
section but different absolute values. In these cases, renormalization based on the most 
reliable absolute cross sections was performed before fitting. 

A few experiments were done with natural Pb targets. In order to compare these data 
with those for pure 208 Pb targets, a search was performed for experimentas in which (p,f) 
cross sections were measured for Pb isotopes [ |3~6| , [3~7| , [3~9"[ ^0[ ]. Then auxiliary approxi- 
mations were determined for the 206 Pb/ 208 Pb and 207 Pb/ 208 Pb proton-induced fission cross 
section ratios (see Fig. 5a, b). The following expression was suggested: 
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Fig. 5. The experimental data [36|, 37, 39[ 5^] and the approximations of the cross section ratios: 
206 Pb(p,f)/ 208 Pb(p,f) (a), 207 Pb(p,f)/ 208 Pb(p,f) (b), na *Pb(p,f)/ 208 Pb(p,f) (c), and na *Pb(n,f)/ 208 Pb(n,f) (d). 
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where r(E) is the cross section ratio, E is the incident particle energy, r , b, E and p are 
fitting parameters found by the least-square method and given in Table III. x 2 l v refers to 
the standard \ 2 P er degree of freedom 

TABLE III 

The parameters of the 206 pb/ 208 Pb and 207 pb/ 208 Pb proton-induced fission cross section 

ratio approximations 



Cross section ratio 


To 


b 


E 


P 


X 2 lv 


206 Pb/ 208 Pb 


1.26 


0.975 


12.8 


1.69 


4.89 


207 Pb/ 208 Pb 


1.15 


16.7 


3.10 


3.74 


3.40 



Combining these approximations with the natural Pb isotopic composition, we deter- 
mined an approximation for the na *pb/ 208 Pb (p,f) cross section ratio (see Fig. 5c). This 
ratio decreases with incident proton energy and becomes nearly constant at an energy of a 
few hundred MeV. 

There are no experiments that include neutron-induced fission cross section measure- 
ments for the relevant separated Pb isotopes within the same experiment. Thus the na *pb/ 208 Pb 
(n,f) cross section ratio was calculated with the use of the nai Pb data of Staples et al. [ |2T| ] 
and the 208 Pb data of Smirnov et al. [ |2"2"| , |2"3| ] (see Fig. 5d). The ratio is constant within the 
experimental uncertainties and is equal to 1.43 ± 0.13 in the 45-160 MeV energy region 
where the data overlap. Assuming that this ratio weakly depends on energy in the 160- 
500 MeV region, which is true for (p,f) reactions discussed above, a renormalization was 
performed for the na *Pb data of Reut [Q, Goldanskiy [|3j, and Staples 0]. 

Figs. 6-9 show the experimental data on the nucleon-induced fission cross sections for 
208 Pb and 209 Bi together with the corresponding approximations described by the following 
expressions: 



a f {E) = exp(c + c x E + c 2 E 2 ) , (2) 
(T f {E) = exp[c + c 1 \nE + c 2 (\nE) 2 ], (3) 
a f (E) = co{l-exp[-ci(E-c 2 )]} , (4) 

where 07 is the fission cross section, E is the incident particle energy, and c , c\ and c 2 are 
fitting parameters found by the least-square method and given in Table IV. 
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cross section. 



14 



102- 



10- 



_Q 

E, 
o 

-4—' 

o 

<D 
(f) 

<Ji 

<D 

O 
i_ 

o 

c 
o 

"w 



1 - 



10-1 - 



10-2- 



10-3- 



10-4 



208pb(p } f) 



♦ 
♦ 



♦ 
♦ 




l 



♦ 



20 



EXPERIMENTAL DATA 

GRADE 1 
■ Shigaev1973 

GRADE 2 

♦ Khodai-Joopari 1966 

• Okolovich 1974 

GRADE 3 

* Kon'shin 1965 
( nat Pb, re-normalized) 

□ Perelygin 1972 



APPROXIMATIONS 



ENDF-B/VI/HE File, 1990 
Prokofyev 1 996 



100 



500 



Projectile Energy (MeV) 



Fig. 8. The experimental data [||, |37|, |jj f|| and approximations §|, [54J] for the 208 Pb(p,f) 
cross section. 



15 



10- 



.2 1 - 



10-1 



208 Pb(n,f) | 



f 



ft 



J* 



,1 



EXPERIMENTAL DATA 

GRADE 1 
O Smirnov 1996, 2°spb 

GRADE 2 
□ Staples, private comm. 
( nat Pb, re-normalized) 

GRADE 3 
A Reut 1 950 

( nat Pb, re-normalized) 
V Goldanskiy 1955 

( nat Pb, re-normalized) 



APPROXIMATIONS 



- ENDF/B-VI/HE File, 1990 
the present work 



20 



100 



Projectile Energy (MeV) 



500 



Fig. 9. The experimental data [ pl| , p3| ], [|21[]-[[23|] and approximations Q54| ] for the Pb(n 
cross section. 



16 



TABLE IV 

The parameters of the nucleon-induced fission cross section approximations for 208 Pb and 

209 Bi 



E (MeV) 


Expression 


c 


Cl 


c 2 


209 Bi(p,f 


) 


37.4-70 


(2) 


-5.150 


0.2203 


-0.001290 


70-3000 


(4) 


182.7 


0.01264 


43.60 


209 Bi(n,f 


) 


20-73 


(3) 


-108.7 


50 


-5.6 


73-1000 


(4) 


100 


0.006 


45 


208 Pb(p,f) 


31.1-70 


(2) 


-16.11 


0.5310 


-0.003760 


70-23000 


(4) 


98.75 


0.01012 


55.05 


208 Pb(n,f) 


30.5-54.5 


(2) 


-4.455 


0.06642 


3.753*10" 4 


54.5-116.8 


(2) 


-4.151 


0.1026 


-3.920*10~ 4 


116.8-462.4 


(4) 


80.71 


3.190*10" 3 


66.06 



The approximations of Fukahori and Pearlstein [ |54| ] included in the ENDF/B-VI/High 
Energy files are also shown in Figs. 6-9. They are seriously discrepant with respect to 
the experimental data below about 40-50 MeV for all four reactions. In the 50-500 MeV 
region there are discrepancies of up to a factor of 2 between Fukahori's (n,f) cross section 
approximations and those given in the present work, which take into account the recent 
experimental data [|2~T|]-[[2"4"1]. The discrepancy for the (p,f) cross section approximations 
is less (no more than 35%). We consider the approximations given in the present work as 
more reliable because they utilize a more complete experimental database and are obtained 
as the result of the consistent procedure of experimental data ranking and selection. 



V. MODELING FISSION CROSS SECTIONS 

The nucleon-induced fission cross section calculations for 208 Pb and 209 Bi nuclei at 
45-500 MeV performed in the present work with the code CEM95 show that there is no 
universal input parameter set that allows us to describe the fission excitation functions in 
the whole energy region considered. In the lower energy part (up to 160 MeV) the best 
agreement with the experimental data is reached with the parameter set presented in Tables 



V and VI. Hereinafter the parameter names are given as in the CEM95 code manual Q20|]. 

Fig. 10 illustrates the results of the calculations with the given parameters in compari- 
son with the experimental data approximations discussed in Sect. IV. 

The calculations reproduce the general trends of the experimental data, although the 
model systematically underestimates the cross sections in the middle of the energy range 
under study (at about 70-100 MeV) and overestimates them at the ends (i.e., above 100- 
150 MeV as well as below 50-60 MeV). An exception is the 208 Pb(n,f) reaction, for which 
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the calculated results lie under the experimental data for neutron energies below 70 MeV. 
This discrepancy may be evidence that the data approximation for the 208 Pb(n,f) cross 
section overestimates the actual cross section at low energies, due to an undetermined sys- 
tematic error in Ref. [|2"T|]. This approximation was deduced in Sect. IV by renormalization 
of the data of Staples [ [ZTj ] assuming that the na *Pb(n,f)/ 208 Pb(n,f) cross section ratio does 
not depend on neutron energy (see Fig. 5d). It is possible that this ratio rises with neutron 
energy decrease below 50-70 MeV, similarly to the na 'Pb(p,f)/ 208 Pb(p,f) cross section ra- 
tio (see Fig. 5c), but the lack of 208 Pb(n,f) cross section data does not allow us to make a 
quantitative estimate. 



TABLE V 

Optimal CEM95 input parameters for nucleon-induced fission cross section calculations 
for 208 Pb and 209 Bi targets in the 45-160 MeV region 



Parameter name and 
recommended value 


Parameter description 


RM=1.5 


The parameter r (fm) in Dostrovsky's formula Q55p 
for inverse reaction cross sections 


WAM (see Table VI) 


The ratio of the level density parameters in the fission 
and neutron emission channels, df/a n 


IFAM=9 


The choice of level density parameter systematics a(Z, N, E*)\ 
the 3rd Iljinov et al. systematics []56|] 


IB=6 


The choice of macroscopic fission barrier model Bf(Z, N): 
Krappe, Nix, and Sierk [57| 


ISH=2 


The choice of shell correction systematics 5W gs (Z, N) for 
the fission barrier: Truran, Cameron, and Hilf [|5E|] 


IBE=1 


The choice of fission barrier dependence on excitation energy 
-B/(.E*):Barashenkov, Gereghi, Iljinov, and Toneev [|27|1 


ISHA=2 


The choice of shell correction systematics 5W gs (Z, N) for the 
level density parameter: Truran, Cameron, and Hilf [|58[] 


IJSP=0 


The choice of fission barrier dependence on angular 
momentum Bf(L): no dependence 



TABLE VI 

Optimal values for the ratio a//a„ in the CEM95 calculations of nucleon-induced fission 
cross sections for 208 Pb and 209 Bi nuclei in the 45-160 MeV region 



Target nucleus 


Incid* 
P 


jnt particle 
n 


209 Bi 
208 pk 


1.21 
1.19 


1.20 
1.20 



18 




50 100 150 200 50 100 150 200 

Projectile Energy (MeV) 



Fig. 10. Comparison between the experimental data and the calculations of the cross sections 
for the following reactions: 209 Bi(p,f) (a), 209 Bi(n,f) (b), 208 Pb(p,f) (c), and 208 Pb(n,f) (d). The 
solid lines represent the approximations of experimental data given in Sect. IV, the dashed lines 
show the CEM95 predictions with the input parameters shown in Tables V and VI. 
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VI. SENSITIVITY OF FISSION CROSS SECTIONS 



Figs. 11-15 illustrate a sensitivity study of the calculated Bi(p,f) cross sections rel- 
ative to variations of the following CEM95 input parameters: 



• the ratio of the level density parameters in the fission and neutron emission channels, 
a f /a n (Fig. 11), 

• the choice of the level density parameter systematics a(Z, N, E*) (Fig. 12a, b), 

• the choice of the macroscopic fission barrier model B f (Z, N) (Fig. 13a, b), 

• the choice of the fission barrier dependence on excitation energy Bf(E*) (Fig. 14a), 

• the choice of the fission barrier dependence on angular momentum Bf (L) (Fig. 14b), 

• the choice of the shell correction systematics 5W gs (Z, N) for the fission barrier 
(Fig. 15a), 

• the choice of the shell correction systematics SW gs (Z, N) for the level density pa- 
rameter (Fig. 15b). 



In all cases, all other parameters are held fixed as each is varied in turn. The calculational 
uncertainties shown in the figures are due to the statistical character of the Monte Carlo 
simulation. 

As seen from Fig. 11, the calculated cross sections are very sensitive to even small 
variations of the af/a n parameter. The variation of the a//a n value by ± 0.03 from the 
optimal value of 1 .21 (see Table VI), about 2.5%, makes the 209 Bi(p,f) cross section vary by 
a factor of 1.5-2, although the shape of the excitation function does not change appreciably. 

Fig. 12 illustrates the calculations with different level density parameter systematics 
[p6|], [|59[]-[|62[]. A comparative study of these systematics and further useful references 
can be found in Ref . JT^ ] . The sensitivity to the choice of the level density parameter sys- 
tematics is relatively weak. The majority of these systematics give nearly the same shape 
to the fission excitation function (with the other parameters fixed according to Tables V 
and VI). The absolute scale of the excitation function varies from one model to another. 
For any model it is possible to improve the agreement with experiment by adjusting the 
a//an parameter. However, it is not possible to alter the shape of the calculated excitation 
function, which is typically too steep above 100 MeV The exception is Malyshev's sys- 
tematics which predicts too gently sloping an excitation function below 70 MeV. The 
shape closest to experiment corresponds to the choice of the 3rd Iljinov et al. systemat- 
ics []5(J, The same choice was found to be the optimal for the CEM in Ref. [JT3p. 

Fig. 13 illustrates the calculational results with different models of the macroscopic 
fission barriers [p7[], y63[j-|J68[j. A comparison of results obtained with these different 
models for fission barriers to experimental data and further references can be found in 
Ref. [[18]]. Similarly to the previous case, the excitation function shape is too steep above 
100 MeV and depends weakly on the choice of the fission-barrier model. The exception 
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is the choice of the constant barrier (curve 8 in Fig. 13b) that leads to a too gently sloping 
excitation function shape below 70 MeV. The absolute scale of the excitation function 
varies typically by a factor of about 2-3 depending on the choice of the specific fission 
barrier model, with the other parameters fixed as given in Tables V and VI. Again it is 
possible to improve the agreement with experiment by adjusting the a//a n parameter. 
The best agreement with experiment was reached with the model of Krappe, Nix, and 
Sierk []57|]. A similar conclusion was made in Ref. JTE|]. 

A similar situation takes place for the choice of the energy dependence of the fission 
barriers (see Fig. 14a). The systematics of Barashenkov, Gereghi, Iljinov, and Toneev [ ]27| ] 
is found to be the optimal one (see Ref. [ |TB| ] for details). 

The choice of the angular momentum dependence of the fission barriers, or more ex- 
actly, the mode of calculation of the saddle-point moment of inertia according to J7(i| ] 
or [ [7T| ] (see details in Ref. [JTBl]) has also a rather weak influence on our present results 
(Fig. 14b). This is due to only a small angular momentum (of only a few K) being im- 
parted to the nucleus by incident nucleons in the energy region under study. Slightly better 
agreement with experiment is reached with no angular momentum dependence of the fis- 
sion barriers. 

Fig. 15 shows the sensitivity of our results to the choice of the shell correction sys- 
tematics [ |58| , |65| , |72| ] for the fission barriers (Fig. 15a) and for the level density parameter 
(Fig. 15b). The shell corrections by Cameron [ f72| ] and Truran, Cameron, and Hilf [ |58| ] 
give similar results, although application of the latter option to the fission barriers (curve 2 
in Fig. 15a) gives somewhat better description of the experiment above 100 MeV. This can 
be explained as follows. The shell corrections from Refs. [|72p and [ |5B[ ] differ significantly 
only for nuclei far from stability. For these nuclei, as was pointed out in Ref. JT8]], the 
corrections Q5B[ ] give a better agreement between calculated and experimentally measured 
fission barriers. Since the fissioning nuclei become more and more neutron-deficient with 
increase of incident particle energy (see Fig. 2), it is natural that Truran, Cameron, and 
Hilf's shell corrections [^HJ] provide a better description of the fission cross sections at 
higher energies. 

Results quite similar to the ones presented in Figs. 1 1-15 were obtained from sensitiv- 
ity studies for other reactions as well. In all cases the model predicts too sharp an increase 
of the excitation functions for energies above 100 MeV. The higher the projectile energy, 
the larger are the discrepancies between the experimental data and the calculations (see 
Fig. 16). 
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Fig. 11. The results of the sensitivity study of the 209 Bi(p,f) cross section calculation with 
CEM95 to the choice of the ratio of the level density parameters in the fission and neutron emission 
channels. The bold line is the approximation based on the experimental data. The thin lines with 
numbers represent the calculations with different a,f/a n . The other model parameters are fixed 
according to Table V. 



22 



100 - 



10- 




± 100- 



10- 



I 

40 



— approximation, 
Prokofyev 1 996 



CEM95 calculations 
with different level density 
parameter systematics: 



a=0.1A 
Malyshev's 
1st Ignatyuk's 
2nd Ignatyuk's 
1st Cherepanov-i 



T 



T 



— approximation, 
Prokofyev 1 996 




CEM95 calculations 
with different level density 
parameter systematics: 

6 - 1st Cherepanov-lljinov's 

7 - 1st lljinov-Mebel's 
2nd lljinov-Mebel's 

9 - 3rd lljinov-Mebel's 

10 - 4th Iliinov-Mebel's 



\ 1 I 1 I 1 I 

100 120 140 160 



60 80 

Projectile Energy (MeV) 



Fig. 12. The results of the sensitivity study of the 209 Bi(p,f) cross section calculation with 
CEM95 to the choice of the level density parameter systematics. The bold line is the approxima- 
tion based on the experimental data. The thin lines with numbers represent the calculations with 
different level density parameter systematics [pq], [59]-[|^]. The other model parameters are fixed 
according to Tables V and VI. 
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Fig. 13. The results of the sensitivity study of the 209 Bi(p,f) cross section calculation with 
CEM95 to the choice of the model for the macroscopic part of fission barrier. The bold line is the 
approximation based on the experimental data. The thin lines with numbers represent the calcula- 
tions with different macroscopic fission barrier models [57], [63]— [58]. The other model parameters 
are fixed according to Tables V and VI. 
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Fig. 14. The results of the sensitivity study of the 209 Bi(p,f) cross section calculation with 
CEM95 to the choice of the systematics for the dependence of fission barrier on excitation en- 
ergy (a) and angular momentum (b). The bold line is the approximation based on the experimental 
data. The thin lines with numbers represent the calculations with different Bf(E*) [27, |6^] and 
Bf(L) [ f70| , [71]] systematics (see [18] for details). The other model parameters are fixed according 
to Tables V and VI. 
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Fig. 15. The results of the sensitivity study of the 209 Bi(p,f) cross section calculation with 
CEM95 to the choice of the shell correction systematics for fission barrier (a) and level density 
parameter (b). The bold line is the approximation based on the experimental data. The thin lines 
with numbers represent the calculations with different shell correction systematics [58, 65, 72] The 
other model parameters are fixed according to Tables V and VI. 
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VII. FURTHER DEVELOPMENT OF THE MODEL 



Since CEM95 allows us to describe well the characteristics of nuclear reactions that 
do not involve fission, like double differential spectra of secondaries (see [Q, H ^]) as well 
as spallation product yields (see [|TT1]-[|T4|]), it is natural to search for the reasons for the 
discrepancies in the modeling of fission. Two different possibilities are: 

1) taking into account dynamical effects in fission, which reflect the connection be- 
tween single-particle and collective nuclear degrees of freedom (see, for example, [f73|]). 



The diffusive character of nuclear motion towards and over the saddle point due to nu- 
clear viscosity leads to a decrease of the value of the fission width. This effect may rise 
with the excitation energy, i.e., in just the direction needed for a better description of the 
experimental fission cross sections; 

2) modification of the calculation of the level density at the saddle point. This possi- 
bility is discussed below. 

As discussed previously, CEM95 incorporates several different level density parameter 
systematics. Most of them utilize the following formula for the excitation energy depen- 
dence of the level density parameter, which was suggested first by Ignatyuk et al. []6l|]: 



a{Z,A,E*) = a{A) 



(5) 



where a is the asymptotic value of the level density parameter at high excitation ener- 
gies 

a(A) = aA + (3A 2/3 B S . (6) 

Here, A, Z, and E* are the mass, charge, and the excitation energy of the nucleus, 
5W(Z, A) is the shell correction, A is the pairing correction, a, (3, and 7 are phenomeno- 
logical constants, and B s is the surface area of the nucleus in units of the surface for the 
sphere of equal volume. For nuclei with equilibrium deformation B s « 1, while for the 
saddle point B s > 1. Formula (5) reflects the effect of the strong correlation between the 
single-particle state density and the shell correction magnitude for low excitation energies, 
and the fade out of the shell effects on the level density for high excitations [EUL FT]] . 



In the original version of CEM95, the level density parameter systematics based on 
formulae (5) and (6) are applied for all decay channels of an excited nucleus except for 
the fission channel. In the latter case, the level density parameter at the saddle point aj is 
calculated using an analogous parameter for the neutron emission channel, a n , and a con- 
stant ratio, af/a n , which serves as a fitting parameter of the model. Thus the shell-effect 
influence on the level density in the neutron emission channel is automatically conveyed 
to the level density at the saddle point. On the other hand, we expect that shell corrections 
at the saddle point for nuclei in this mass range should bear no relation to those at the 
ground state, due to the large saddle-point deformation, and the consequent different mi- 
croscopic level structure near the Fermi surface . In fact, the shell corrections at the saddle 
point should be of much smaller magnitude than those at the ground state for nuclei in 
the neighborhood of 208 Pb, due to the greatly reduced symmetry of the saddle-point shape 
compared to the spherical ground state (see also Refs. [ J73] , |73p). 
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In order to estimate the importance of taking into account this effect, we perform cal- 
culations with the parameter a/ being energy-independent, which is equivalent to the dis- 
appearance of the shell-effect influence on the level density at the saddle point. The mag- 
nitude of the parameter aj was calculated using formula (6) with the same values of the 
coefficients a and (3 that were utilized for the other decay channels. The parameter B s 
was adjusted to provide the best agreement with the experimental fission cross section 
approximations (see Sect. IV). All the other model parameters were the same as for the 
calculations described in Sect. V (see Table V). The optimal values of the parameter B s 
are shown in Table VII. 

TABLE VII 

The optimal values of the parameter B s for the fission cross section calculations with the 

modified version of CEM95 



Target nucleus 


Incid( 
P 


mt particle 
n 


209 Bi 
208 pb 


1.18 
1.15 


1.12 
1.12 



The values for the parameter B s are physically reasonable because they are larger than 
1, which reflects the larger deformation of the fissioning nucleus at the saddle point in 
comparison with the equilibrium state and they are not far from the value B s ps 2 1 ' 3 , 
which is expected for the saddle-point configuration of the preactinides [[76]]. However, 
we do not yet understand why different values are needed for proton- and neutron-induced 
fission, and why the values of B s in the neutron channel are significantly smaller than those 
corresponding to the known deformation of saddle-point shapes in this mass region. 

The ratio a//a^, where is the asymptotic (large E*) level density parameter in the 
neutron channel, is fixed by the fitting of B s , and is not independently fitted as was done 
previously, is calculated using formula (6) and the values a = 0.072 and [3 = 0.257, 
which correspond to the 3rd Iljinov et al. systematics [|56l]. The ratios appear to depend 
only weakly on the nuclear mass. For the reactions under study their values (given in 
Fig. 16) are in the range 1.045-1.070, which are consistent with the compilation of 5//a^ 
values published by Ignatyuk et. al. [f76|]. 

The results of calculations with the modified CEM95 are shown in Fig. 16. We see a 
much better description of the experimental fission cross sections in comparison with the 
original version for nucleon energies of 100-500 MeV On the other hand, the calculation 
systematically overestimates to a slight extent the fission cross sections below 100 MeV. 
We emphasize that other related improvements need to be made before a final assessment 
of the model's value and predictive capability can be made. For example, an expected 
excitation- energy dependence of the ground- state shell correction would change the aver- 
age height of the calculated fission barriers as the incident energy is increased. However, 
the current modifications do indicate the crucial importance of properly incorporating ap- 
propriate level densities and motivate the search for a consistent model of barriers, ground- 
state masses, and level densities which may improve the predictive power of the CEM. 
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Fig. 16. Comparison of the experimental data on the 2U9 Bi(p,f) (a), 2uy Bi(n,f) (b), 2US Pb(p,f) (c), 
and 208 Pb(n,f) (d) cross sections and the calculations with the original and modified versions of 
CEM95. The solid lines represent the approximation of the experimental data (see Sec. IV). The 
dotted lines show the original CEM95 results with the parameters given in Tables V and VI. The 
dot-dashed lines show the results after the modification of the fission channel in CEM95. 
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VIII. SUMMARY 



We have performed a systematic analysis of the nucleon-induced fission cross sections 
for 209 Bi and 208 Pb nuclei in the 45-500 MeV energy range using the CEM95 code, which 
has incorporated several choices for the level density parameters, nuclear masses, shell 
and pairing corrections, saddle-point moments of inertia, and for the macroscopic and 
microscopic fission barriers. We have performed a detailed analysis of the dependence of 
calculated fission cross sections and of distributions of the residual nuclei after the cascade 
stage of reactions and of fissioning nuclei at the compound-nucleus stage of reactions on all 
these choices. We find that the distributions of residual and fissioning nuclei with respect 
to their mass, charge, and excitation energy are not very sensitive to the specific model 
choices used in our calculations. Also, our results agree satisfactorily with the scanty data 
available in the literature for the average share of the linear momentum transferred to the 
fissioning nuclei. 

We present here analytical approximations determined by a critical analysis of all avail- 
able experimental data, and compare them to calculations using CEM95 with and without 
modifications. Our analysis shows that by choosing the appropriate value for the ratio of 
the level density parameters in the fission and neutron emission channels a,f/a n , it is pos- 
sible to describe satisfactorily fission cross sections of all reactions under investigation in 
a limited range of energy. But it is impossible to get a satisfactory agreement in the whole 
range of energy using a fixed set of CEM95 parameters and of nuclear models used in the 
calculations. This means that the version of the CEM as realized in the code CEM95 does 
not allow us to predict arbitrary fission cross sections in a large range of incident energies. 
We perform in the present work a further development of the model, where we have mod- 
ified the calculation of the level density parameter of nuclei at the saddle point, by use of 
simple physical considerations. This change allows us to get much better agreement with 
the experimental fission cross sections in the energy region from 100 to 500 MeV, but fur- 
ther improvement of the CEM95 is still needed at energies below about 100 MeV. We plan 
to continue our work by extending the range of incident energies and of nuclear targets 
investigated, striving for a model capable of predicting fission cross sections for arbitrary 
targets in a wide range of incident energies. 
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